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Abstract We derive thermal QCD sum rules for the correlation function of two vector currents in the rho-meson 
channel. It takes into account the leading non-perturbative corrections from the additional operators, which appear 
due to the breakdown of Lorentz invariance at finite temperature. The mixing of the new operators has a drastic 
effect on their coefficients. The thermal average of all the operators can be related to that of the quark condensate 
Qs^ • and the energy density. The sum rules then yield the temperature dependence of the parameters of the p-meson, 
namely its mass and coupling to the vector current. Our result is that these parameters are practically independent 
of temperature at least up to a temperature of 125 MeV. 
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I The QCD sum rules Q], proposed about two decades ago, prove remarkably successful in addressing the non- 
perturbative problems of hadron phenomenology. In this approach one considers the product of two local operators, 
like the currents of the QCD theory. A sum rule is obtained by equating the dispersion relation for its vacuum matrix 
element at sufficiently large space-like momenta to the corresponding Wilson operator product expansion The 
higher dimension operators present in the expansion give rise to the non-perturbative, power corrections. The idea 
of extending these sum rules to finite temperature by replacing the vacuum matrix element with the thermal average 
naturally suggests itself. 

O 1^ The original work establishing the thermal QCD sum rules is that of Bochkarev and Shaposhnikov They 
O , recognised the importance of the (low energy) continuum states in the spectral function to account for the effects of 
' the medium. On the basis of these sum rules they discussed the temperature dependence of the resonance parameters 
^ , and the existence of phase transition. However, there arise additional operators in the Wilson expansion at finite 
■ temperature [Q, which were not correctly incorporated in their sum rules: In effect, they calculated these new operator 
rS I contributions perturbatively, which cannot be justified, particularly at low temperature. 

' The additional operators arise because of the breakdown of Lorentz invariance at finite temperature by the choice 
of the thermal rest frame, where matter is at rest at a definite temperature The residual 0(3) invariance naturally 
brings in additional operators. The expected behaviour of the thermal averages of these Lorentz non-invariant (new) 
operators is somewhat opposite to those of the Lorentz invariant (old) ones: While the old operators start with 
non-zero values at zero temperature and decrease in magnitude with the rise of temperature, the new ones, on the 
other hand, are zero at zero temperature but grow rapidly with temperature. The importance of including these new 
operators in the thermal sum rules, particularly at not too low a temperature, is thus clear. 

Although a number of works on thermal QCD sum rules exist by now, these are flawed with respect to the new 
operators: Either some of these are missed |^ or their mixing, which changes their coefficients drastically, is not taken 
into account ||^. In a recent work Q we applied a simple, configuration space method Q, to evaluate the Wilson 
coefficients of these new operators (up to dimension four) which appear at finite temperature in the short distance 
expansion of the product of two quark bilinear operators. Here we make use of this result to write and evaluate the 
thermal QCD sum rules, incorporating correctly the contributions from all the dimension four operators. 

We consider the correlation function of the time ordered (T) product of two vector currents in the p-meson channel. 
The use of the T-product, rather than the retarded (or advanced) product, is a little complicated in writing down the 
spectral representation but has the advantage in perturbative calculations, for which we can apply the conventional 
formalism. Throughout this work we shall employ the real time formulation of the thermal field theory , which 
requires in general not only the physical fields but also the accompanying 'ghost' fields. Since, however, we work to 
lowest order in perturbation expansion, ghost fields do not show up. 
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It is convenient to write difference sum rules by subtracting the vacuum sum rules from their finite temperature 
counterparts. For such sum rules the absorptive parts are expected to be saturated with the p-meson pole and 
the TTTT-continuum ||^, whose contributions we derive here for completeness. The thermal averages of the different 
operators reduce essentially to that of the quark condensate and the energy densities of quarks and gluons fl^. Chiral 
perturbation theory has been used to calculate the temperature dependence of the quark condensate |14|. The 
difficulty with the energy densities is that while one of them is the total energy density, which can be obtained from a 
knowledge of the hadronic spectrum at low temperature, the other is an unphysical combination of the two densities. 
We need an additional input to relate this latter combination to the total energy density. 

Once the power corrections are known, the difference sum rules give the temperature dependence of the p-meson 
parameters, namely its mass and its coupling with the vector current. With our simple saturation scheme, the sum 
rules can be used up to a temperature of about 125 MeV. The numerical evaluation shows that these parameters have 
negligible dependence on temperature. 

In sec. II we write the kinematic decomposition for the thermal correlation function of two vector currents, derive 
the Landau representation for the time ordered product and state the results of operator product expansion to 
derive finally the form of the thermal QCD sum rules. In sec. Ill we obtain the contributions of p and tttt intermediate 
states to the spectral function. In sec. IV we collect the information on thermal average of the operators present in 
the sum rules and evaluate the difference sum rules. In sec. V we discuss how to extend the sum rules to higher 
temperature. In the Appendix we give the details of the evaluation of a limit stated earlier Q . 



II. SUM RULES 



We derive the QCD sum rules for the thermal average of the time ordered (T) product of two currents, 

V'i<i) - * / d'xe^"-- {T {V;{x)V^{o))) . (2.1) 

Here V^(a;) is the vector current (in the p meson channel) in QCD theory, 

V^{x) = q{x)-i^{T^/2)q{x), (2.2) 

q(x) being the field of the u and d quark doublet and t"" the Pauli matrices. The thermal average of an operator O 
is denoted by (O), 

(O) = Tre-'^"0/Tre-^", 

where H is the QCD Hamiltonian, (3 is the inverse of the temperature T (coinciding with the time ordering symbol !) 
and Tr denotes the trace over any complete set of states. 



A. Kinematics 



At finite temperature Lorentz invariance is broken by the choice of a preferred frame of reference, viz, the matter 
rest frame where temperature is defined. But the book-keeping with indices becomes simpler if we restore it formally 
by introducing the four- velocity of the matter l^. (In the matter rest frame = (1, 0, 0, 0).) The time and space 
components of g^, are then raised to the Lorentz scalars, lo = u ■ q and q = \J uS^ — . We shall, however, return to 
the matter rest frame while doing all actual computations. 

In such a Lorentz invariant framework, there are two independent conserved kinematic covariants fl^ , which we 
choose as 

P -_ Mil^^~ ~ n - ~ 

where U/j, — u^j, — Luq^/q^ . The invariant decomposition of T^^ is then given by 

V'il) = S'''{Q^,,Ti + P^.Tt), (2.3) 
where the invariant amplitudes Tj.t are functions of q^ and lu. 
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Notice that the kinematic covariants P^i^ and are free from singularities at = (and also at = ). This 
is convenient at finite temperature as there are dynamical singularities extending up to = 0. With this choice of 
kinematic covariants, the dynamical singularities reside only in the invariant amplitudes. 

To extract the invariant amplitudes from T^^, it is converient to form the scalars Hi. 2, 

Then the invariant amplitudes are given by 

q^ 2 \ q 

In the limit q = 0, the amplitudes T; and Tj are related. To see this we write the spatial components of 7^;^ as 



Ti = ^Ha, Tt^--[Ui + ) . (2.4) 



'r,f{q)^5-%,Tt-q,q,{Tt~qlTi)], 
where qi is the ith component of the unit vector along q. As ^ 0, there cannot be any dependence on getting 

Ttiqo, \q\ = 0) = qlTiiqo, \^ = 0). (2.5) 



B. Spectral representation 

Let us obtain the spectral representation for the correlation function in qo at fixed \q\. First, evaluate the trace 
over a complete set of eigenstates of four-momentum, when it becomes a sum over forward amplitudes weighted by 
the corresponding Boltzmann factors. Then insert the same set of complete states between the currents to extract its 
X- dependence which is then integrated out. Introducing a 5- function in qo, the result can be written as 

27r yq'^-qa-ie q'^ - qo + le J 

where the expression for with the double sum over states may be converted back to the form, 

M^liq) = [ d'xe^'^^V^{x)V^{o)). (2.7) 



Using the double sum representation, it is easy to show that 

M^liq,, k1) = e-^*M;^(go, k1). (2.8) 

The opposite sign of it in the two terms in eqn.(2.6) is typical of T- products. As a result the imaginary part of 
is given by the sum, \{M^^ + M^°), while the principal value integral contains the difference, ^{M^^ — M^°). 
But the two are related by (2.8), 

M^^t - Mil = + ^^'D tanh(/3go/2). (2.9) 

We thus get the Landau representation for the T product at finite temperature 

r;.'(9o, k1) = ^ImT^'Aq^^ \^) + P j^^ dg^,%^«OlMl, (2.10) 

% - 90 

where, for brevity, we write 

iV;^(g) = S'^'N^M = n-'lmT;^{q) tanh(/3go/2) (2.11) 

Expressing q'^ = q^, + (90 ^ 'Zo)w^t, we recover the representation for the invariant amplitudes Ti^t- Further, using 
the symmetry properties, ImTi t{—q^) = ImTi t{qf_i) and going over to imaginary values of go (^o ~ ^Qcn Qo 
for which A'^jy vanishes, these become, 

"'0 9o + 

It may actually require subtractions, but it docs not affect the Borcl transformed sum rules we shall write below 
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C. Operator product expansion 



The contributions of operators of dimension four to Ti^t are obtained from the short distance expansion and improved 
upon by the renormahsation group equation in a previous work Including the unit operator (the perturbative 
contribution), T; is given for large Euclidean momenta by (Q^ — —q^), 

^ In(^) + (m^u) + + I^{(e> + aiQ^)iWef/3 - 9^) ■ (2.13) 

Tt is also given by the same expression, except for an overall factor of —Q^ and a factor of (1 — 2(f /Q'^) multiplying 
the term with 0's. Here m is the degenerate mass of the u and the d quarks and ^{qq) — {uu) ~ {dd). The operator 
is quadratic in the gauge field strength GJ^^,, (a = 1, • • • , 8), = (Q:s/7r)G°^G°'^'', with Us — .g^/47r, where g is the 
gauge coupling constant. Along with these two operators already contributing to zero temperature, we now have the 
linear combinations of two new ones, O^^'^ = u^^Q-^^'^u'^ , where Qf^'^ are the energy-momentum tensors for quarks of a 
single flavour and of gluons. The corresponding components for the total tensor is — UfQ^ + O^, where is the 
number of effective quark flavours. The logarithmic dependence of a{Q^) arise due to the anomalous dimension d 
of the operator (16e-''/3 - 0^), 

where 1 GeV) is the scale at which all renormalisations are carried out. 

Note the mixing of the operators O-'^ and in (2.13) under the renormahsation group, which is well-known in 
the context of deep inelastic scattering iQ. In the operator product expansion of two quark currents, the Wilson 
coefficients of Q-^ and are, to leading order, of zeroth and first order in as, arising from Born and one-loop 
graphs respectively. However, due to operator mixing, the coefficients are drastically changed in that both 6 and 
(160^^/3 — ) have coefficients with leading terms of zeroth order in a^. In (2.13) we retain only these leading terms. 

D. Sum rules 

We now equate the spectral representation and the result of operator product expansion for the amplitudes T/ and 
Tt at sufficiently high Qq. Taking Borel transform we arrive at the thermal QCD sum rules For Ti we get 

where (O) is the non-perturbative contribution of higher dimension operators, 

(o) = Muu) + ^ + Y^q^{(Q) + «(M')(i6eV3 - e«)}. (2.15) 

and a similar one for Tf. Each is a two parameter sum rule, dependent not only on T but also on \q\. 

In the thermal rest frame the thermal average of the new operators are energy densities, which increase rapidly 
with temperature. Earlier works on thermal QCD sum rules were flawed, as these contributions were not properly 
included. 



III. ABSORPTIVE PARTS 

We work below the critical temperature, where hadrons constitute the physical spectrum. As with the vacuum sum 
rules, the dominant contribution to the spectral function is given by the p- meson. We also calculate the contribution 
of the non-resonant TTTr-continuum. The question of other significant contributions will be discussed in Sec. IV. 
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A. p-pole 



The coupling of the vector current to the p-meson is given by 

<Q\V^\p^>=5'^''Fpmpe^, (3.1) 

where is the polarisation vector of p of mass mp. The experimental value of Fp as measured by the electronic decay 
rate of p° [g is Fp = 153.5 MeV. 

A simple way to calculate the absorptive part of the p-pole diagram is to note the field-current identity of 
with the rho- meson field pj^, 

V^{x)^mpFpp1{x), (3.2) 
which reproduces (3.1). Then the p-meson contribution is given essentially by its thermal propagator, 

rpl\q)^tmlF'p J d^xe^'^^p1Xx)pl{o)) 

= ^<5'^X^,^(^-ry,. + ^)A^,(g), (3.3) 

where l^1^{q) is the 11- component of a scalar field propagator with mass nip, 

A?i(<7) = ^y— + 2^n{uj,)6{q^ ~ ml), (3.4) 



and n{ujq) is the Bose distribution function, niuq) — {e^^i — 1) ^ , ujg = yep + m^. We then get 

Loop corrections at finite temperature will make rrip and Fp temperature dependent, modifying them, in general, 
differently in the longitudinal and transverse amplitudes. These modifications may be obtained by calculating the 
appropriate loop graphs. Here we shall find them by evaluating our sum rules. 



B. TTTT-continuum 



The TTTT-contribution to the amplitudes describes the interaction of the current with the particles in the medium, 
which are predominantly pions. Chiral perturbation theory gives the contribution of the pion field (ff-ix) to the 
vector current, which, to lowest order, is 

V^{x) = e^'''<p\x)dp<P^{x) 
Then the pions contribute to the correlation function as 

T^'il) - ^5'^' j (01(2^ - 9)^(2^ - 9).Aii(fc)Aii(fc - (3.6) 

where AJ^ is again the 11- component of the scalar propagator (3.4) but with mass m^r. Its imaginary part can be 
obtained by the cutting rules at finite temperature [|l^. Here we obtain it directly for this simple amplitude. We 
express as 

A^i(fc) - {1 + n{uk)}D{k) + n{uok)D*{k), (3.7) 

where D[k) is the zero temperature propagator, D{k) — ijik^ — rn^ + ie), and carry out the ko- integration pO[| . The 
imaginary part may then be read off as 

ImT^^iq) - S-''{Lp,{q) + Lp,{-q)}, (3.8) 
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where 

I d^k {2k~q)^{2k~q), 



(27r)3 Aujiuj2 
[{(1 + ni)(l + 712) + nin2}(5(qo - ^^i - ^^2) 
+ {(1 + ni)n2 + (1 + n2)ni}(5(qo - t^i + ^^2)]- (3.9) 



Here ni = n[uji), n2 = f^(w2) with ui = y k'^ + , uj2 = y {k — q)^ + rri^. The time component of k^ in the tensor 
structure is understood to be given by fco =1^1. 

With the help of the 5-functions we can rerwite (3.9) to get N^^ defined by (2.11) for go > 0, 



d^k (2fc-g)^(2fc-g)^ 
(27r)3 4wiW2 

[(1 +ni + n2)S{qo ~ uji - L02) + (n-2 - "■i)^('7o ~ ^2)]. (3.10) 



In this form the factors involving the density distributions can be interpreted in terms of pion absorption from and 
emission into the medium Since the first and the second 5- functions in (3.10) contribute to time-like {q^ > Arri^) 
and space-like (q^ < 0) regions respectively, we write it as 

/d^k 
^^^(2fc-,),(2^-,). 

[(1 + ni + n2)e{q^ - Ami) + («2 - ni)ei-q^)]S{iqo ~ LUi)^ - cui). (3.11) 

Thus in addition to the usual cut, 4m^ + l?!^ < 9o < 00, the amplitude at finite temperature acquires a new short 
cut, < g2 < kT ^,§. 

The angular integration is carried out using the 6- function, when the constraint |cos0-^| < 1 results in a 9- 
hmction, 9[—q^{uJi — a;+)(wi — w_)], where 

1 



w± = ^ 1^^")' " v^l - 4m2/g2. (3.12) 



We thus get 



1 r°° dioi 



■(^ZK - q)f_,[ZK - q)^[L + 711+ n2)a[ 

{2k-q),,{2k^qUn2-ni)e{-q^). (3.13) 



2|g1 (2^)2 

It is now simple to extract the absorptive parts of the invariant amplitudes by using (2.4). Changing the variable 
uji to X given by uji — ^{qo + \q\x), we get 



for g2 > 4^2 (-3 -^4^ 



with 



and 



°° / 0^2 



A^ry 64^2 (72(^2 _ 

{n{Mx-qo)/2)~n{Mx + qo)/2)}, for < q (3.I6) 

The superscripts (±) on denote time- like and space-like q^^ respectively, where they are non- vanishing. The first 
term on the right of (3.14) arising from the unity in the factor (1 -f tii -I- 712) in (3.13) is the zero temperature 
contribution of the tttt state. Evaluated here in a non-covariant way it, of course, agrees with the covariant evaluation 
of the Feynman amplitude (3.6) with ^^^{k) replaced by the vacuum propagator D{k) [l^. 
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C. Explicit sum rules 



Let us now write explicitly the sum rule (2.14) for T;. Saturating Ni with the above contributions it becomes, 

Jiml + \q\'^ Jo ) 

As the temperature goes to zero, the two terms in bracket go to zero and the thermal average of the operators on the 
right become the vacuum expectation values, recovering the familiar vacuum sum rule. The integral on the left is the 
non-resonant 27r contribution and is small compared to the resonance contribution. 

As \q\ 0, the sum rule (3.17) simplifies considerably. The limit for the second integral in bracket is given in ref. 
Here we derive it in the Appendix. The sum rules for T; and Tt then become, 

J.2(y)^-mJ(T)/M^ + /o(Af2) + It{M') = £ + (3.18) 



and 



where 



Io{M^) = 


487r2 , 


POO 

Jiml 




487r2 , 


poo 

/ dsse-'^^'\\ 

J Ami 


/y(M2) ^ 


247r2 , 


/>oo 

/ ds{e-''^^'\^ + v{i - v^)/2}n{^s/2), 

J 4m2 


Jj,(M^) = 


247r2 , 


/ dsse-''^'\^n{y/:s/2), 

J 4m J 



(3.20) 

with V = v{s) = Y^l — 4m2/s. Observe that the sum rules (3.18-19) are not independent, in agreement with the 
relation (2.5): the second one is obtained by differentiating the first with respect to l/A/2. 



IV. EVALUATION OF SUM RULES 



In the above sum rules we have approximated the absorptive parts of the amplitudes by the p-meson pole and the 
2TT continuum, while we retain only the contributions of the unit operator (the perturbative result) and of all the 
operators of dimension four in the operator product expansion. To check this saturation scheme, let us compare the 
zero temperature limit of our sum rules with the corresponding vacuum sum rules [0. The latter include, in addition, 
the rather large contribution from the high energy continuum beyond 1.5 GeV, as indicated by the experimental data, 
as well as the contribution of a quark operator of dimension six, which is also large because of its origin in Born 
(rather than loop) diagram. Since we do not include any of these contributions, wc cannot expect the sum rules as 
written above to be well saturated. 

Rather than incorporate these contributions, we isolate the thermal effects by considering the difference sum rules, 
obtained by subtracting out the vacuum sum rules from the corresponding finite temperature ones. Then the con- 
tribution to the absorptive parts beyond 1.5 GeV, being temperature independent, cancels out in the difference. 
Thus it is the temperature dependent contributions of the p-meson and of the 2tt continuum, which should domi- 
nate the absorptive parts of these sum rules. Also the dimension six quark operator, Og say, contributes an amount 
^ (Oe) — (OlOelO), which is insignificant in the immediate neighbourhood of T = and as the temperature rises, this 
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contribution is overwhelmed by that of the two-gluon and other Lorentz non- invariant operators, as our estimates 
below for these operators show. 

The difference sum rules for the two invariant amplitudes allow us to calculate the temperature dependence of the 
/9-meson parameters, 

Amp(T) EE mp(r) - nip 



\^--i^^'y-^^^M^\^^)\^ (4.1) 




\ m?, I 



with 



(O) = rnluu) + ^ + -^^^m + «(M2)(16eV3 - 6^)}, (4.3) 

where the bar on the operators indicates subtraction of their vacuum expectation values. Here we insert the experi- 
mental values for mp and Fp, as these are well reproduced by the vacuum sum rules. 

We now collect information on the operator contributions. The vacuum expectation value of the chiral condensate 
(0|gg|0) is known from the PCAC relation of Gell-Mann, Oakes and Renner, 

2to(0|um|0) = -F^ml, 

where m = ^{rriu + mj) and the pion decay constant i^^r, defined by, 

{0\A^^\7:\q)) =zq^d^''F^, 

has the value = 93.3 MeV. Taking m = 7 MeV (22|, we get (0|uu|0) = -{225MeV)^. The vacuum expectation 
value of the two-gluon operator, as determined from the QCD sum rules 0, is (0|G^|0) — (3301/6^)*. 

The operator is related to the trace of the energy momentum tensor Q^^, by the trace anomaly. Normalising it 
to zero vacuum expectation value and taking thermal average, it gives 



(G2) = (G^) - (0|G2|0) 

--^(<©m)-E"^^(^)- (4-4) 

The trace at finite temperature is given by (6p = (6) — 3P, where (Q) is the energy density and P the pressure. 

The temperature dependence of both (uu) and (0) have been calculated in chiral perturbation theory [|4j . Cor- 
rections due to nonzero quark masses as well as the contributions of the massive states {K,r],p, ■ ■ •) have also been 
incorporated. However, as the authors point out, the validity of the perturbation theory and the use of dilute gas ap- 
proximation to calculate the contribution of the massive states restrict these results to within a temperature of about 
150MeV. Thus the critical temperature Tc is, strictly speaking, beyond the range of validity of their calculation. 
Since, however, the order parameter (uu) falls rapidly at the upper end of this range, one has only to extrapolate it 
a httle further to get Tc = 190MeV. 

Besides the total energy density (9), there also occurs the thermal average, (lee-'^/S-ef). The last one cannot be 
calculated without further input, at least in the hadronic phase. Now both naive counting of the degrees of freedom 
and empirical study of the pion structure function Q suggest the quark fraction of the energy density to be about 
half of the total. So we assume 71/(6-'') = x/(0),with Xf = -5, whence (16e-''/3 - O^) = {(^ + l)xf - l}(e). 
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FIG. 1. Shift in the rho-meson mass as a function of temperature for M'^ = 1 GeV'^ and = 4 GeV^. 

As with the zero temperature sum rules, the results are expected to be stable in a region in Af^, which is neither 
too high to make the continuum contribution large relative to the resonance contribution nor too low to emphasize 
the neglected power corrections of higher order. Since the high energy continuum contribution gets cancelled in 
the difference sum rules, the region of may be extended somewhat at the upper end. Figs. 1 and 2 show the 
evaluation for equal to IGeV^ and AGeV^. The results for Airip and AFp are rather stable for temperatures 
up to about 125 MeV. At higher temperatures the results, particularly for ArUp appear unstable. Closer observation 
reveals here a large cancellation between the 27r contribution and the leading power correction we have retained. Thus 
the non-leading power correction become important here, whose inclusion may restore the stability in to higher 
temperatures. 




25 50 75 100 125 150 

T (MeV) 

FIG. 2. Shift in the coupling of rho-meson with the vector current as a function of temperature for = 1 GeV^ and 
=4 GeV^ 



V. DISCUSSIONS 

In this work we have written the thermal QCD sum rules for the two point correlation function of vector currents in 
the p-channel, including the leading power correction due to all the operators of dimension four. Because of the loss of 
Lorentz invariance at finite temperature, two new (Lorentz noninvariant) operators creep up in the Wilson expansion, 
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in addition to the two old (Lorentz invariant) ones, already existing in the vacuum sum rules. Thus compared to the 
two numbers, (0|g(7|0) and (0|G^|0) in the vacuum sum rules, we have now four temperature dependent quantities, 
{qq), (G^), (G) and {160^/3 — 8^) in the thermal sum rules. All these thermal averages can be evaluated from chiral 
perturbation theory, supplemented by contributions from the massive states. An extra input is needed only for the 
last unphysical operator. The sum rules can then be used to determine the temperature dependence of the p- meson 
parameters. Our evaluation shows that the mass of the p-meson and its coupling with the vector current remain 
practically unaffected by the rise of temperatur e up to at least 125 MeV. The absence of the shift in the mass appears 
to agree with one set of results obtained in ref [p3| . 

Unfortunately the sum rules, as they stand, cannot be extended up to the critical temperature. One reason 
contributing to this restriction has to do with the evaluation of the thermal average of the operators, as we already 
discussed in the last section. What restricts it further is their instability under a change of for temperatures 
above 12bMeV. Even in this restricted temperature range, the numerical evaluation shows that the new operators 
are significant. In fact, for temperatures above say 70 MeV, the new contributions overwhelm the old ones in the 
difference sum rules. This situation necessitates reanalysis of all earlier results based on thermal QCD sum rules, 
where the new operators are not properly taken into account. 

We now discuss a possible way to extend the sum rules to higher temperatures. In the vacuum sum rules the 
operators of dimension four (and higher) provide corrections to the leading (perturbative) result. But in the difference 
sum rules it is these corrections which become the leading contribution. One thus expects that by including nonleading 
contributions from higher dimension operators along with those of dimension four already included, the sum rules 
would be stable against variations in over a wider range of temperature. This inclusion is all the more necessary 
for sum rules like the one for ArUp^ where the leading operator contribution cancels largely with that of TTTr-continuum. 

The higher dimension operators will, of course, complicate the evaluation of the sum rules in that we have to know 
the temperature dependence of their thermal averages. Also there is more proliferation of operators than what is 
generally thought. The procedure in the literature |^ of including dimension six quark operators and excluding the 
Lorentz noninvariant gluon operators |2j], because of the smallness of their coefficients by a factor of as(Af^), is 
not justified. For, as we have seen, the quark and the gluon operators mix under a change of scale, so that after 
renormalisation group inprovement, both the coefficients are of the same order in as(M^). 

A way to proceed here is to write the entire set of sum rules by considering two-point functions of not only the 
vector quark bilinear but also the others, like the scalar, tensor, etc. All of these sum rules receive contribution from a 
few resonances and the operators from the same set. The simultaneous evaluation of all these sum rules is expected to 
provide a self-consistent check on the thermal average of the operators. Further, using quark bilinears of appropriate 
chiralities, one can get sum rules without any of the gluon operators |^5|| . These sum rules should prove easier to 
evaluate and would also check the saturation scheme in a simpler context. 
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Appendix 



Here we calculate the limit of the second integral in bracket in (3.17), 

1 flo]'^ 1"°° 

^^32^ Jo ^9oe-^'°"''''^^*' y dxx'{nmx~qo)/2)-nmx + qo)/2)} (A.l) 

as |g| ^ 0. It is convenient to change the integration variables qo and x to A and u respectively by go = M'll ^-^^ 
\q\x = u, getting 



^^shl! rf^'e"^''"'^''^'/'"' duu'^{nm^ ~ 9o)/2) - n((|g> + qo)/2)}, I = ^ \g]^ + ^ 
As |g| — > 0, it becomes 



IGtt^ Jq J -2^s= d,u 

1 



487r2 



dsv{-i-v^)n{y/^/2) (A.2) 

4m?. 
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To get the second line we integrate by parts over u and interchange the order of integration in the remaining double 
integral 

The hmit of the corresponding integral for Tt is zero, because of the presence of in the integrand. 
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